Topological and topological-electronic correlations in amorphous silicon 



o 
o 

(N 
Oh' 

o 



CO 



c 

o 
o 



> 
in 

o 
d 

o 



X 



Yue PanQ Mingliang ZhangQ and D. A. Drabokfj] 

Department of Physics and Astronomy, Ohio University, Athens, OH 45701 
(Dated: February 5, 2008) 

In this paper, we study several structural models of amorphous silicon, and discuss structural and 
electronic features common to all. We note spatial correlations between short bonds, and similar 
correlations between long bonds. Such effects persist under a first principles relaxation of the system 
and at finite temperature. Next we explore the nature of the band tail states and find the states to 
possess a filamentary structure. We detail correlations between local geometry and the band tails. 

PACS numbers: 61.43.Dq, 71. 23. An, 71.23.Cq 



I. INTRODUCTION 

The nature of three-dimensional space-filling disor- 
dered networks is salient to a variety of topics in scientific 
research. Despite work in recent decades, a number of 
puzzles remain both about the structure and connectiv- 
ity of such networks. Also, important quantities derived 
from such structures (such as electronic or vibrational 
states) are quite incompletely understood, though their 
characteristics are ultimately derivative from the under- 
lying structure. In this paper, we systematically consider 
a collection of (largely) tetrahedral continuous random 
networks of amorphous Si, and and despite the fact that 
several of these are quite "standard", we detect and re- 
port novel features of their connectivity and electronic 
structure. 

The character of the tail states in disordered semicon- 
ductors is a problem of importance with a history dating 
back at least to the sixties. It has long been understood 
that band tail states in amorphous semiconductors arise 
from strains in the network sufficient to push states past 
the band edges into the gap. Such tails decay exponen- 
tially into the gap: this feature, nearly universal to dis- 
ordered semiconductors, is called Urbach tailing^ Clever 
photoemission experiments allow the separate observa- 
tion of the valence and conduction tails and even the 
temperature dependence of the tailing^. Several theories 
have been presented to explain the Urbach tails, most 
recently by Cohen and coworkers 

Beside band tailing and the detailed structure of the 
networks we study, our work is related to the general 
problem of the topology of space-filling networks. For 
example, work on the cuprates^^ provides clear, albeit 
circumstantial evidence for the existence of percolating 
filamentary networks, and is seen in STM measurements 7 
(presumably two-D projections of the three-D filamen- 
tary network). Filamentary paths are also believed to ex- 
ist in fast ion conducting glasses with larger than average 
free volume^ regions providing pathways for ion trans- 
port (see for example Fig. 3 in Ref£). Farther afield, 
results akin to ours (the percolation of short bonds) can 
be inferred in protein binding in a yeasA 

This paper has two immediate aims. First, considera- 
tion of several models reveals certain persistent geometri- 



cal features, notably strong self-correlations between con- 
figurations involving short bonds and analogous correla- 
tions involving long bonds. To the extent that these cor- 
relations exist in various structural models and survive 
ah initio thermal molecular dynamics (MD) simulations, 
it seems worthwhile to point these out. The second aim 
is to discuss general features of the band tail states in a- 
Si. Resonant cluster proliferation of tail and gap states is 
observed^, and filamentary connections between local- 
ization centers. In agreement with earlier report a 11 ' 12 , 
we observe a robust tendency for valence (conduction) 
tail states to be associated with short (long) bonds. Un- 
like our earlier report^, we did detect analogous features 
for the bond angle distribution, where we see that smaller 
(larger) bond angles are correlated with valence (conduc- 
tion) tails. 

These general observations about the tail states in a- 
Si are relevant to calculations of transport, since it is 
presumably just these states through which carriers hop, 
as we discuss in detail elsewhere^. 



II. CONSTRAINTS ON THEORY 



Model characteristics and nomenclature 



Even in good quality unhydrogenated material, mid- 
gap defects are fairly rare, at most a few sites per 
thousand-^. At a practical level this means that ideal 
models of tractable size (up to several thousand atoms) 
can contain zero to a few defects producing mid-gap 
states. The Wooten-Weaire- Winer (WWW) 15 method 
is the "gold standard" for forming such realistic mod- 
els, and several of the models we discuss depend upon 
WWW at least as a starting point. For completeness, 
and in particular to ascertain the extent to which our 
observations might tend to be WWW artifacts, we have 
repeated the calculations for other models as well. We 
summarize some salient features of the models in Table 
I. 
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TABLE 1. Models and basic information. NOA 
stands for "Number of atoms" ; De% stands for "de- 
fects percentage"; A0/Acos# stands for the width 
of bond-angle 9 /cos 9 distribution; Ab stands for the 
width of bond-length distribution; WWW stands for 
Wooten-Weaire- Winer modeling scheme^; ART stands 
for activation-relaxation technique^; PC stands for 
"paracrystalline" ; RMC stands for Reverse Monte-Carlo 
modeling scheme^. 

Mi, M2, M3, Ms and Mg all derived originally from 
WWW 15 modeling scheme. Mi, M2 and Mg were devel- 
oped by B. R. Djordjevic, M. F. Thorpe, F. Wooter*i£; 
M 3 was developed by Feldman and M 8 by N. Mousseau. 
M4 was made via activation-relaxation technique^. M5 
and M 6 are Reverse-Monte-Carlo^ models. Model Mi 
to M4 and Mg to Mg are continuous random network 
(CRN) models. Beside CRN models, we include model 
M5, a paracrystalline model (CRN with crystalline inclu- 
sion). In particular for M5, 211 crystalline silicon atoms 
are embedded. 

Excepting M5, all the models that we have studied 
here have normally distributed bond lengths and cosines 
of bond angles. For the large WWW models the sim- 
ilarity to the normal distribution is quite striking 21 . 
Here we report the bond length and bond angle dis- 
tributions for these models in Fig. [TJ Model Mi and 
M4 show nearly perfect Gaussian fits, while model M7, 
(a 500-atom model), is fairly Gaussian, with statistical 
noise. The other models exhibit similar traits. Dong and 
DraboldSi have attempted to correlate these normally- 
distributed quantities to Urbach tailing. 

Model Mi, M 5 , M 8 and M 9 are defect free. In the 
other models, there are coordination defects, ranging 
from 0.2% to 12.5%, including both dangling bonds and 
five-fold "floating" bonds. The pair correlation functions 
of these models are all similar, accurately reflecting the 
experimental function. As usual with a-Si, this is a re- 
minder that the information in the pair-correlation func- 
tion alone is very insufficient to specify coordinates. 



B. Local inter-bond correlations 

Our inspection of the topology of the models reveals 
that there is a tendency for short bonds to be linked to 
other short bonds and long bonds to be linked to other 
long bonds. In particular, this leads to 1) a tendency to 
spatial separation of the shortest and the longest bonds 
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FIG. 1: Bond length distribution for model Mi, M 4 and M 7 
(left column); Cosine bond angle distribution for model Mi, 
M4 and M7 (right column). 

in the system; 2) a tendency to build up chain and ring 
like structures among extreme bonds. Selected fractions 
of shortest and longest bonds of model Mi are extracted 
from the models and displayed in Fig. f2J The correlation 
is clearly visible, especially in Fig. [2};-f. Interpenetrating 
filamentary structures appear to percolate through space 
for a sufficiently large cutoff in the fractions of bonds dis- 
played. For comparison, the 4% shortest and 4% longest 
bonds of model M2 - M7 are shown in Fig. [3] To some 
extent, the self-correlations of short and long bonds are 
visible in every model shown, suggesting that such bond 
correlations are not an artifact of a particular model or 
modeling scheme. 

To explore our observation more quantitatively, we 
computed bond-bond correlation functions for the short- 
est and longest bonds: 

y Nl N 2 

«i = l n 2 fni 

where N±, N 2 are the number of short or long bonds 
as needed; n\ and n 2 count over these particular bonds; 
r nin2 is the distance between the bond centers of bond 
rii and bond n 2 ; V is the volume of the unit cell. The 4% 
shortest and 4% longest bonds are taken for model Mi 
and the corresponding bond-bond correlation functions 
for (1) Ni = 4% shortest, N 2 = 4% longest; (2)N t = 4% 
shortest, N 2 = 4% shortest; (3) Ni = 4% longest, N 2 
— 4% longest are plotted in Fig. The first peaks of 
the correlations of the same-type bonds (grey and dot- 
ted lines) appear to be much stronger than the correla- 
tion between the different- type bonds (solid black line). 
The second and even the third peaks maintain this fea- 
ture, and after 6.0 A, this correlation wanes. This con- 
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FIG. 2: (A) 1%, (B) 2%, (C) 3%, (D) 4%, (E) 5% and (F)8% 
shortest (dark) and longest (light) bonds of model Mi. 
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FIG. 4: Bond-bond correlation function between selected 
bonds of model Mi. Solid black line denotes pair correla- 
tion function between the 4% shortest and the 4% longest 
bonds. Grey line denotes pair correlation function among the 
4% shortest bonds. Dotted line denotes pair correlation func- 
tion among the 4% longest bonds. 
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FIG. 3: 4% shortest (dark) and longest (light) bonds of model 
M 2 , M 3j M 4 , Ms, M 6 and M 7 . 



variations from model to model. Model Mi, M3, M5 and 
Mg show slightly stronger effects than the rest. Model 
M3 and M5 show stronger correlations between the short 
bonds compared to those between the long bonds, while 
model Mg and M7 have relatively stronger correlations 
between the long bonds than those between the short 
bonds. 



C. Persistence of correlation effects with ab initio 
relaxation 

The correlations are found quite consistently in the 
several models we examined, suggesting that the effects 
are not modeling artifacts. We also undertook a sim- 
ple comparison to see how a conjugate gradient (CG) 
relaxation process by SIESTA affects the models. A re- 
laxation was done on model Mi and Mg. While there 
were slight modifications in structure due to the CG pro- 
cess, we didn't observed any significant changes in these 
correlations. 



D. Persistence of correlation at finite temperature 



firms that there is a correlation among the same-type- 
bonds, with a correlation radius around 6.0A for model 
Mi. The very similar correlation functions between the 
4% shortest and the 4% longest bonds of model M2 - M7 
are plotted in Fig. [5] A substantial correlation between 
same-type bonds is revealed for all models, with minor 



To further explore the likelihood that the long and 
short inter bond correlations were artifacts, we also un- 
dertook thermal MD simulations at 300K on two small 
models (M 8 and M 9 ) with the code SIESTA. Nose dy- 
namics were used. With some thermally-induced fluc- 
tuations, the short-short bond correlation and long- long 
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FIG. 5: Bond-bond correlation functions between selected 
bonds of model M 2 , M 3 , M 4 , M 5 , M 6 and M 7 . Solid black 
line denotes pair correlation function between the 4% shortest 
and the 4% longest bonds. Grey line denotes pair correlation 
function among the 4% shortest bonds. Dotted line denotes 
pair correlation function among the 4% longest bonds. 

bond correlation persist at all times, as seen in Fig. [5] for 
model Mg, where the heights of the first peaks of the pair 
correlation functions from Eq. (fTJ) are plotted throughout 
time steps. 4% shortest bonds and 4% longest bonds of 
the system at all time were involved for the calculations. 
The correlations between the same types of bonds are 
fluctuating but generally higher than that between the 
different types. The detailed temperature dependence in 
this correlation, however, requires further study. 



III. ELECTRONIC STRUCTURE 



FIG. 6: Correlation (1) between the 4% shortest bonds and 
the 4% longest bonds; (2) among the 4% shortest bonds; (3) 
among the 4% longest bonds of model Mg throughout time 
steps at 300K. Nose dynamics were used. 
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A. Density of states and Localization 

The electronic and optical properties of the models are 
determined by the coordinates of the atoms. Electronic 
eigenstates and wavefunctions of the models are obtained 
for static lattices. For the larger model M2 and M4 a tight 
binding mode l 22 i 23 is applied. For other models, the ab 
initio code SIESTAS! is used. Results are very similar 
for the empirical and ab initio calculations. 

Densities of states near E/ of the seven models are 
given in FigJTl Spatial localization is reported using the 
Inverse Participation Ratio (IPR) in Fig. [5J Model Mi 
and M5 are defect free, with no states in their optical 
gaps. Model M 2 has 0.2% defects and one gap state. 
Model M3 and M4 both have 4% defects, and a small 
peak of states in the gap. M 6 has 12.5% defects and 
M7 has 12% defects and both have many states spread 
through the gap. Only Mi, M2 and M5 conform closely 



FIG. 7: Densities of states of models. Dashed lines indicate 
optical gaps. Fermi levels are all shifted to 0.0 eV. 



to ordinary expectations for an a-Si DOS. M3, M4, Mq 
and M7 exhibit too many defects for good a-Si materials 
and are studied here mainly for comparison purpose. 



B. Topology of band tail states 

1. Filamentary structure 

The most- localized states of model Mi, displays a 1-D 
filament connectivity. For state #1024(-0.34eV, a valence 
state) and State #1025(0. 34eV, a conduction state) the 
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FIG. 8: Spatial localization obtained from IPRs (Inverse Par- 
ticipation Ratio) of models. Dashed lines indicate optical 
gaps. Fermi levels are all shifted to 0.0 eV. For M2 and M4 
only a small part of IPRs near the Fermi level is shown. 



most electron-probable atoms are selected and shown in 
Fig. [HI Interestingly, in both case, the atoms locate 
closely (localization) and happen to form filament (chain) 
like structures. In addition, the angles along the chain 
of E#1024 are all small angles, and the bonds are all 
short bonds. The angles along the chain of E#1025 are 
all large, and the bonds all long. 

Localized states of the other six models, M2 to M7, 
were also examined. We picked the most localized state 
for each model. Only the atoms possessing most charge 
are shown in Fig. We see that, for all these mod- 
els, these localized atoms form filament-like structures in 
general. The strength of the effect varies only slightly 
among the models. 



2. Bond length- electron correlations and bond 
angle- electron correlations near Ef 

Fedders, Drabold, Nakhmanson's studyii observed 
that valence tail states preferentially involved short bond 
lengths, whereas conduction tail states tended to be in- 
volved with longer bond lengths. Angle distortion was 
also mentioned but no specific relation to the band tails 
was detected. Related work was undertaken in Rcf. 12 . 
Here, we extend this analysis to our collection of mod- 
els. Where the bond lengths are concerned, the present 
work is consistent with these findings. In addition we 
also clearly detect a correlation between valence tail and 
small bond angle, and a correlation between conduction 
tail and large bond angle, which has not been explicitly 
observed before. 

An electron(Mulliken charge)-weighted mean bond 
length for each state was calculated in Reb^i, and it re- 



FIG. 9: Most localized electronic eigenstates of model Mi, 
M 2 , M 3 , M 4 , Ms, M 6 and M 7 . Approximately 60% - 
75% charge of the state is shown in each case. For model 
Mi, both a valence tail state(#1024) and a conduction tail 
state(#1025) are shown. 



vealed an interesting and general trend in bond lengths. 
Here we computed the mean bond length in a similar 
fashion, and explored into seven of our models for this 
mean value. We then extended the calculations into the 
mean bond angles, using the very same Mulliken charge. 

For a certain eigenstate E, the Mulliken charge q( n ,E) 
is the electron probability on an atom n. Using it as 
a weighting factor, we define the corresponding sym- 
metrized mean bond length as B( E ^ . 



B, 



J2, hm b (n, m )q(n,E)q{ m ,E) 



(2) 



where n, m summations only go over all the possible 
bonds in a unit cell (of N atoms), with n, m as the bond 
end atoms. q( n ,E), as defined, is the electron probabil- 
ity on the atom n. (near the Fermi level) for each 
model is plotted in Fig. [TTJ] B/e) shows an asymmetrical 
split around the gap (Fermi level) for all models. Smaller 
average bond lengths are correlated to the valence tails 
and large average bond lengths correlated to the con- 
duction tails. There are clear difference among the seven 
models, especially on the mid-gap states from the models 
with high defect concentration. For models with fewer 
defect states in the gap, the bond-electron correlations 
are very clear. Models with many mid-gap states show 
mixed behavior, though the asymmetric splitting is al- 
ways detectable. 

In analogy with the discussion for bond lengths, we 
define the electron- weighted mean of bond angles as: 
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FIG. 10: B(E), the electron- weighted average of bond length 
near the Fermi levels for the seven models. Only small frac- 
tions of data are showing for M2 and M4. 



FIG. 11: A(E), the electron-weighted average of bond angle 
near the Fermi levels for the seven models. Only small frac- 
tions of data are showing for M2 and M4. 



where n,m,l are the atoms of the unit cell of N atoms; 
summations of n, m, 1 only count over all the possible 
bond angles, with m as the vertex of each angle. q( n ,E)i 
as defined, is the electron probability on the atom n. At e) 
is plotted in Fig. [TTJ It also shows a clear asymmetri- 
cal splitting around the gap for all plots. Smaller average 
bond angles clearly appear to be correlated to the valence 
tails and larger average bond angles correlated to the con- 
duction tails. There are also differences among the mod- 
els but the asymmetric feature remains. Interestingly, 
these angle-electron correlations show strong resemblance 
to the bond-electron correlations discussed above. It is 
also natural to connect these correlations to the filamen- 
tary structures shown in Fig. [9l For Mi, which is a 
typical example, the structure from state #1024 (a va- 
lence tail state) has both short bonds and small bond 
angles along the chain, while #1025 (a conduction tail 
state) has both long bonds and large bond angles along 
the chain. 

One feature to note is that here only the energies near 
the Fermi level are shown. In our models, A^e) an d Bt e) 
also have extrema in other spectral gaps, whenever the 
states are fairly localized. 



IV. DISCUSSION 

Here we offer some heuristic explanation of the last two 
sections. In a-Si, there are two kinds of random poten- 
tials leading to the localization of carriers: point defect 
(missing atoms or impurities) and topological disorder. 
The former produces very localized states which lie in 
the neighborhood of the Fermi level. These are bound 
states around the point defects, and sit in the middle of 
gap. The scattering suffered by the carriers is so strong 



that the localization length of those states is smaller than 
the average bond length (for example, the lone electron in 
dangling bond). 25 In contrast, topological disorder, the 
deviation of bond lengths and bond angles from their 
crystalline values, produces localized states primarily ly- 
ing in the tails of the valence band and the conduction 
band. The localization length of those tail states may 
extend from one bond length to a range including many 
atoms depending on the structure of the network. 

Band tails exhibit localization because the density of 
states in their energy range is small, and there is less 
opportunity for derealization via mixing, as discussed 
in the "resonant cluster proliferation" modeh 10 i 26 In this 
picture, it is doubtful that a an initially localized energy 
state could persist in the midst of a continuum; it would 
inevitably mix with the reservoir of resonant states, if 
there was any overlap between the reservoir and localized 
state. Similarly, a well-localized state would survive if 
placed in a region of zero density of states. From this 
point of view, the band tails are simply between these 
limits. 

According to the scattering theory of solidS ) 25 ' 27 i 28 the 
distortion of the atomic configuration from crystal mixes 
Bloch states with different band indices and wave vec- 
tors. If these secondary scattering components have large 
enough phase shifts relative to the primary Bloch wave 
outside the distorted region, they will interfere destruc- 
tively, because the phase memory related to the primary 
Bloch wave has been lost. While inside the distorted 
region, secondary scattering components are almost in- 
phase with the primary Bloch wave (phase shift is less 
than 7r), and they will interfere constructively inside this 
region. Thus if a distorted region can produce scattered 
components with large momentum transfer, a localized 
state can be formed by the interference of secondary com- 
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ponents with the primary Bloch wave. 

With this mechanism, it is possible that several lo- 
calized states have peaks in the same distorted regions. 
This was observed in previous work^S on M2. We ex- 
plored other models, M4 and M7, and found similar states 
with overlapping clusters in space. On the other hand, 
one localized state may have several peaks in structurally 
distorted regions which are spatially separated (Fig. [9]). 

The observation that valence tail states are preferen- 
tially localized on short bonds and conduction tail states 
localized on long bonds (Fig. [TU)) can be understood as 
follows. Short bond length has two consequences: (1) it 
will increase the transition integral; the energy of valence 
band states are lowered while the energy of conduction 
states is lifted relative to a hypothetical reference crys- 
tal (diamond for a-Si); (2) Small bond lengths increase 
the charge at the center of the bonds (states near va- 
lence edge). Thus if we make the significant assumption 
that the valence and conduction tails should be asso- 
ciated with either short or long bonds, it is clear that 
the electronic (band-energy) is optimized if the valence 
tails are associated with short bonds, and the long bonds 
would "do the least damage" if put above the Fermi-level, 



or linked with the conduction tail. This completely ne- 
glects other contributions to the total energy, though the 
observation of the effect in density functional calculations 
suggests that the band energy is the key term. 



V. CONCLUSION 

We have discovered significant and reasonably model- 
independent features of the network and electronic struc- 
ture of a-Si. We believe that the results may prove rele- 
vant to a number of other systems as welli^ 
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